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The stability of the Bianchi type I anisotropic brane cosmology is analyzed in this paper. We also 
study the efltect of the brane solution by comparing the models on the 3-brane and the models in 
the conventional Einstein's space. Analysis is presented for two different models: one with a perfect 
fluid and the other one with a dilaton field. It is shown that the anisotropic expansion is smeared 
out dynamically for both theories in the large time limit independent of the models with different 
types of matter. The initial states are, however, dramatically different. A primordial anisotropic 
^SJ ' expansion will grow for the conventional Einstein's theory. On the other hand, it is shown that the 

. initial state is highly isotropic for the brane universe except for a very particular case. Moreover, 

' it is also shown that the Bianchi type I anisotropic cosmology is stable against any anisotropic 

' perturbation for both theories in the large time limit. 
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I. INTRODUCTION 



> 

00 I The observation of the cosmic microwave background (CMB) radiation |2| indicates that our Universe is globally 
OO , isotropic to a very high degree of precision. Therefore, our Universe is usually assumed to be described by the 
Friedmann-Robertson- Walker (FRW) metric in most of the literatures. The origin of the isotropic universe has also 
become an interesting research topic ever since. 

On the other hand, it is known that there is a small large-angle CMB anisotropies, AT/T ~ lO^'"^, under the CMB 
] background 0, |^ . The isotropy of our Universe may have to do with the choice of initial conditions and the stability 
of the evolutionary solutions. Therefore, the present isotropic phase could be a dynamic result of the evolution of our 
Universe, no matter what the initial state started out. In this paper, we will address on this issue in both the standard 
JL Einstein's theory and the newly developed brane world scenario 31. In particular, we will consider the evolution 
^ of an anisotropic cosmology described by the Bianchi type I models gJ-[^ for two different types of matter sources: 
(~| , one is the perfect fluid and the other one is a dilaton field. Our result shows that the anisotropic models considered 
■ all evolve dynamically from the anisotropic Bianchi type I universe into the isotropic FRW space in the large time 
limit. This property reveals that the isotropy of the cosmological principle may be justified and made consistent with 
our current observational data. 

5h I These two different theories give, however, completely different initial anisotropic expansion at the very early stage 
of evolution. Indeed, for the conventional Einstein's theory (CET), the anisotropic expansion tends to be large in the 
very early stage. In another words, the universe has to begin from a highly anisotropic initial expansion and then 
decays to zero as the time increases. On the other hand, in the brane world scenario due to the quadratic 

correction which significantly changes the early time behavior of the Universe. As a result, any initially non- vanishing 
anisotropy parameter tends to vanish in the very early period. There is a characteristic time, ^c, that divides the 
evolution of mean anisotropy parameter into two different stages. The mean anisotropy parameter is increasing when 
t < tc and reaches its maximal value a,i t = tc- After that, mean anisotropy parameter starts to decay. This result 
remains true for both the model with a perfect fluid and the model with a scalar field. And this appears to be a 
general feature independent of the types of matter considered. 

In addition, the stability analysis indicates that all models in two different theories are all stable against 

any anisotropic perturbation. In particular, we also show that the system is stable against the dilaton perturbation 
when the dilaton potential is close to its local minimum. Note that the exact solution is known only when the dilaton 
potential is a constant. This also justifies the stability of the isotropic universe in the large time limit under the effect 
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of a dilaton potential with local minima. 

This paper is organized as follows. In Sec. II we briefly review the brane world formulation and the Bianchi type I 
model. Then we consider the evolution of the anisotropic cosmology and its stability in both theories of the GET and 
the brane world. Analysis is presented for model coupled to a perfect fluid in section III. The model with a dilaton 
field and a constant potential is presented in section IV. Finally, the conclusion and possible implications are drawn 
in Sec. V. 



II. PRELIMINARY OF BRANE WORLD AND ANISOTROPIC COSMOLOGY 



The brane world scenario assumes that our Universe is a four-dimensional space-time, a 3-brane, embedded in the 
5D bulk space-time. All the matter fields and the gauge fields except the graviton are confined on the the 3-brane 
as a prior requirement in order to avoid any violations with the empirical results. Moreover, inspired by the string 
theory/M-theory, the Z2-symmetry with respect to the brane is imposed ||2^. A formal realization of the brane world 
scenario, which recovers the Newton gravity in the linear theory, is the Rundall-Sundrum model ^ in which the 4D 
flat brane(s) is embedded in the 5D anti-de Sitter (AdS) space-time. Later on, a covariant formulation of the effective 
gravitational field equations on the 3-brane has been obtained via a geometric approach by Shiromizu, Maeda and 



Sasaki Il3, 15 . It is shown that the effective four-dimensional gravitational field equations on the brane take the the 



following form 

where G^^ and T^^^ are the Einstein and energy-momentum tensors. Sf^^ is a quadratic contribution of T^^ defined as 

^^"^ ^ ul^^^" ~ i^^""^""" ^ 2i^'"' (3r"^ra/3 - T^) . (2) 

The effective 4D parameters, e.g., the cosmological constant A and gravitational coupling ^4, are determined by the 
5D bulk cosmological constant A5, the 5D gravitational coupling and the tension of the brane A via the following 
relations 

A = fci(f + ^j, kl = ^. (3) 

Here g^^ is the metric tensor on the brane. In addition, the quantity E"^^ is a pure bulk effect defined by the bulk 
Weyl tensor 

^From the Eq. (|l|), it is easy to realize that the brane world scenario is different from the GET by two parts: 
(a) the matter fields contribute local "quadratic" energy- momentum correction via the tensor S*^^, and (b) there 
are "nonlocal" effects from the free field in the bulk, transmitted via the projection of the bulk Weyl tensor E^^. 
Therefore, the GET can be treated as a limit of the brane theory by taking E^^, ^ and fcs ^ with properly 
adjusted values of the constants fc4 and A. 

In addition to the generalized Einstein equations (|l|), the energy- momentum tensor also satisfies the conservation 
law VpT^"^ = 0. Therefore, there is a constraint on the tensor iJ^j^, ^^Ef^" = k^S/^^S'^'', due to the Bianchi identity 
on the brane. Here, the operator V is the covariant derivative with respect to the metric g^^^, on the brane. One should 
point out here that the field equations on the brane, namely the generalized Einstein equations, the conservation of 
energy-momentum and the constraint on E^i, are, in general, not a closed system in the 4D brane since the quantity 
E^i, is five-dimensional. It can only be evaluated by solving the field equations in the bulk. In this paper, we will, 
however, only consider the quadratic effect on the brane world in the anisotropic background. Therefore, we will set 
E^i, = which is equivalent to embedding the 3-brane in the pure AdS bulk space-time. 

Another important subject in this paper is the anisotropic cosmology described by Bianchi type I metric. The line 
element of the Bianchi type I space, an anisotropic generalization of the flat FRW geometry, can be written as 

ds'^ = -dt^ + a\{t)dx^ + al{t)dy^ + al{t)dz^ , (4) 

with ai{t), i = 1, 2, 3 the expansion factors on each different spatial directions. For later convenience, we will introduce 
the following variables 

3 

V = JJ^Ci, volume scale factor, (5) 

i=l 
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Hi = — , 1 = 1,2,3, directional Hubble factors, (6) 

ai 

1 -A V 

H = - 2, Hi — -Tj, mean Hubble factor. (7) 
In addition, we will also introduce two basic physical observational quantities in cosmology: 



{Hi — H) , 

A = y ——75 , mean anisotropy parameter, (8) 

^ — ' 3il^ 

1=1 

q = '^(H^^) — 1, deceleration parameter. (9) 

Note that ^ = for an isotropic expansion. Moreover, the sign of the deceleration parameter indicates how the 
Universe expands. Indeed, a positive sign corresponds to "standard" decelerating models whereas a negative sign 
indicates an accelerating expansion. 

III. ANISOTROPIC UNIVERSE WITH A PERFECT FLUID 

In this section, we will consider the case that the matter energy- momentum tensor, T^jy, is a perfect fluid whose 
components are given by 

T% = -P, T\ = T\ = T% = p. (10) 

Here the energy density p and the pressure p of the cosmological fluid obey a linear barotropic equation of state of 
the form p = (7 — l)p with 7 a constant in the range 1 < 7 < 2. 

A. Conventional Einstein's Theory 

For the GET, the dynamics of the space-time is determined by the Einstein equations and the energy-momentum 
conservation law 

G^. = -Ag^, -I- kjT^,, V^.T'^'' = 0, (11) 

with Gpjy the Einstein tensor, A the cosmological constant and fc4 the gravitational coupling A:| = SttG. For the 

Bianchi type I cosmology with a perfect fluid the tt- and w-components of Einstein equations ^ and the conservation 
law become 

3H + j2H? = A-^-^kjp, (12) 

1=1 

:^^(^if.) = A-^fc2^, z = l,2,3, (13) 

p + 3-/Hp = 0. (14) 



First of all, the Eq. (14) can be easily solved to obtain the time evolution law of the energy density of the fluid 

p = poV'"^ , pq = constant > 0. (15) 
Next, by summing the Eqs. (|l3|) with respect to the index i, we can obtain the relation 



^ Actually, the equations and ( |l^ are the time and spatial components of the following version of the Einstein equation: Rij,u 
^g^v + k'^{T^v — \Tg^^), with R^i, the Ricci tensor and T the trace of T^^. 
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FIG. 1: Mean anisotropy parameter of the Bianchi type I conventional Einstein's theory with a perfect fiuid: 7 = 2 (solid 
curve), 7 = 1.5 (dotted curve) and 7 = 1 (dashed curve). The normalization of the parameters is chosen as 3A = 1, 3fc|po = 2, 
and C = 1. 



Subtracting this result with the Eqs. (|T^), we can show that 

V ' 



z= 1,2,3, 



(17) 



with Ki, i = 1, 2, 3 the constants of integration satisfying the consistency condition Y^^=i ~ 0- Moreover, by using 
the evolution law of the matter energy density ( [l5| ) , the basic equation ( p^ describing the dynamics of the anisotropic 
Universe can be written as 



(18) 



with the general solution 



where function G{V) is defined by 



G{V) = 3A1/2 _^ 3fc2^py2-7 ^ (J 



(19) 



(20) 



Here C are constants of integration. Based on this result, the other variables can be calculated straightforwardly and 
the answers are 



H = 



GiV) 



1/2 



3V ' 



Ki / V-^G{V)-^^^dV 



, z = l,2,3. 



A = 3K^G{V)-\ 

IS-fklpoV'^"'' + 12G 
~ 4G(F) 



(21) 

(22) 
(23) 
(24) 



Here the ogi, K^, i = 1,2,3 are constants of integration and = Y^^'f^iKf. In addition, the arbitrary integration 
constants Ki and G must satisfy the consistency condition 



- -G. 
3 



(25) 



Although the general solution can only be expressed in the parametric form via the volume scale factor V, the 
behaviors of the physical variables, namely the mean anisotropy parameter A and the deceleration parameter q, can 
be plotted with respect to the time as shown in the Fig.|l| and Fig.^j. 

Moreover, from the the above parametric expressions, we can still analyze the physics in two different interesting 
limits. First of all, let us consider the large t (or equivalently large V) limit. Since the value of parameter 7 is in 
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FIG. 2: Deceleration parameter of the Bianchi type I conventional Einstein's theory with a perfect fiuid: 7 = 2 (solid curve), 
7 = 1.5 (dotted curve) and 7 = 1 (dashed curve). The normaUzation of the parameters is chosen as 3A — l,3fc|po ~ 2, and 
C= 1. 

the range 1 < 7 < 2, the asymptotic value of G{V) will approach G{V) 3AV^ when the cosmological constant is 
non- vanishing (A > 0). This leads the volume scale parameter of our Universe to expand exponentially in the large 
time limit, i.e. V oc exp[-\/3At]. Therefore, in this limit, the mean anisotropy parameter decays to zero exponentially, 
A cx exp[— 2\/3At] 0, and the deceleration parameter becomes negative q (x exp[— 7V3At] — 1 < 0. Hence the 
Universe can be isotropized dynamically and undergoes an accelerated expansion in the large time limit due to the 
presence of a positive cosmological constant A. In fact, the value of the cosmological constant can only change the 
expanding rate of the Universe but will not affect its effect of isotropization in general. 

In order to illustrate this point, let us consider the model with a vanishing cosmological constant A = 0. Therefore, 
we have G{V) cx V^^''' in the large time limit. This leads to the resuhs V cx t^^f, A cx ^^-4/7^ and q 87/2 - 1 > 0. 
Hence, the Universe can still be isotropized dynamically except for the case when 7 = 2. Note that the evolution will, 
however, be deaccelerated in the case 7 = 2. 

Next, let us focuss on the earlier stages of the above exact solutions. For simplicity, we will assume C = in this 
case. One can show that G{V) ^ V'^~^ . Therefore, one can solve for V cx i^/'' from the Eq. (|2l|). Hence one has 
A cx t^""*/"*", and 87/2— 1 >0 from the rest of the field equations. Therefore, at the earlier period stage of the 
Universe, the evolution is decelerating even a cosmological constant is present. This is also shown numerically in the 
Fig.|l| and Figj^. Moreover, the mean anisotropy parameter A is always non-vanishing independent of the values of 7. 
This means that the early Universe is always anisotropic. Therefore, for the Universe with perfect fluid matter in the 
GET, the initial state is always anisotropic and this primordial anisotropy is smeared away as a consequence of the 
evolution of the Universe. 



B. Brane Cosmology 

We will focuss on the brane effect for the model with a perfect fluid in this subsection. The same constraint for 
the perfect fluid conservation law (^4|) still holds on the brane. In addition, the gravitational field equations and the 
conservation law on the brane take the form, in terms of the variables 

8ij + ^ i/,^ = A - ^ fc> - ^ kip\ (26) 

ll(VH,) = A-2-lklp-2^ky, . = 1,2,8, (27) 
p + 3jHp = 0. (28) 

By comparing the above equations with the ones we considered in the previous section for the GET, it is easy to 
realize that the difference is the quadratic effect due to the energy density p. Moreover, the brane cosmology will 
reduce to GET if we take the limit ^ and adjust the value of ^4 accordingly. We will perform the stability 
analysis, in the next section, for the brane universe. The GET can be reproduced by imposing the limit k^ — > 0. 

The general solution of the above system was obtained in the exact parametric form in [l^ by the same approach 
used in the previous section. Instead of expressing results as functions of time, we are able to present all variables. 
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including time, in terms of volume scale factor, V, with V > 0. For instance, the time variable can be expressed as 

t-to^ I F{V)-'/^dV, (29) 



where F{V) is defined as 



F{V) = 3AV^ + 3kjpoV^-'> + -klplV^-^^ + C, 



with po and C the constants of integration. The other variables are 

F(y)i/2 



H = 

0.1 

A 

q 



aoiV^/^ exp 



Ki j V-^F{V)-^l'^dV 

m^F{v)~^, 

4i^(F) 



1,2,3, 



-1, 



(30) 

(31) 

(32) 

(33) 
(34) 



where ooi, Ki, i = 1,2,3 are constants of integration and — ^^^iKf. In addition, the arbitrary integration 
constants Ki and C must satisfy the same consistency condition K^ = 2C/3. 

One can immediately show that the effect of the energy density quadratic term becomes significant at the high energy 
epoch, or in another words, at the early stages of the Universe by looking at the Eq. (pO|). Indeed, F{V) oc y^"^''' 
at the limit < when V is extremely small. As a result, the F diverges as f — > 0. Therefore, the mean anisotropy 
parameter A ^ at the early universe. On the other hand, in the large time limit, the properties of the Universe 
should be more or less the same as the case we have discussed for the GET in the previous section by looking at the 
same Eq. (^0|) . Hence the early evolution of the anisotropic Bianchi type I brane Universe is dramatically changed 
due to the presence of the brane correction terms proportional to the square of the energy density. The time variation 
of the mean anisotropy parameter of the Bianchi type I space-time is presented, for different values of 7, in Figj^. 
From the Fig.|[ it is clear that at high energy density the evolution of the brane Universe always starts out from an 
isotropic state with A —> 0. The mean anisotropy parameter increases and reaches a maximum value after a finite time 
interval tc- One can show that, when t > tc, the mean anisotropy parameter is a monotonically decreasing function 
approaching zero in the large time limit. This behavior is in sharp contrast to the usual evolution in the GET, as 
shown in Figj^, in which the Universe has to kick off from state of maximum anisotropy due to the constraint from 
the field equation. 

In addition, the early time evolution of the brane universe is normally not in an inflationary phase. On the 
other hand, the brane Universe always ends up in an inflationary phase in the large time limit in the presence of a 
nonvanishing cosmological constant. These are generic features of the brane Universe due to the constraint of the 
field equations on the brane cosmology. Indeed, a more detailed information can be extracted from the Eq. (|2^) in 
the limit t 0, ov equivalently, the case with a small V. For simplicity, one will take C = again. Indeed, we can 
show that V oc i^^'^ as i — * 0. Hence, the expansion of the early universe is of the form of power law expansion. In 
addition, in the early stages of evolution of the brane Universe the mean anisotropy parameter varies as A oc t^~'^^'^ 
approaching zero as i ~ 0. Moreover, the deceleration parameter is given by g = 37 1 which is always positive for 
all possible values of 7 for the case C = 0. 



C. Stability Analysis 

The general perturbations for the FRW background with perfect fiuid can be found in Ref. jl^. The same 
consideration is, however, more complicate for the anisotropic background. For the primary effect, we will only 
consider the scalar mode and neglect the vector and tensor modes [|18| , |l9| . The metric perturbation considered here 
is 

ai ^ QBi + Sai = aBi{l + Sbi), (35) 
and the perturbations with respect to the perfect fluid considered in this paper is 

Pb +Sp, P ^ Pb + {1 - 1)15/0. (36) 
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FIG. 3: Mean anisotropy parameter of the Bianchi type I brane Universe with a perfect fiuid: 7 = 2 (solid curve), 7 — 1.5 
(dotted curve) and 7=1 (dashed curve). The normalization of the parameters is chosen as 3A — l,3fc|po = 2, fc|po ~ 4, and 
C= 1. 
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FIG. 4: Deceleration parameter of Bianchi type I brane universe with a perfect fluid: 7 = 2 (solid curve), 7 — 1.5 (dotted 
curve) and 7 = 1 (dashed curve). The normalization of the parameters is chosen as 3A — 1, 3fc|po ~ 2, fc|po — 4, and C = 1. 



Here the variables with subscript B are the exact solutions presented in the previous sections. For technical conve- 
nience, we will use variables Sbi instead of Sui in our analysis. Therefore, the perturbations of the following quantities 
can be shown to be 

H,^HB; + 5k, H ^ Hb + ^Y.^^^^ Y.^? ^Y.^B^ + ^Y.^B^Sb,, V -> Vb + ^ 56. • (37) 

i i i i i 

As discussed in the previous sections, the results in the GET are special cases for the brane universe, in the limit 
^5 = 0. Therefore, we can perform our stability analysis only for the brane universe first. The conclusion for the GET 
can be easily extracted from the brane world case. 

The perturbation equations for various quantities, can be obtained by substituting the perturbations ( p6| , |37| ) into 
the field equations (^ |2^, ^). Leading terms will reproduce the field equations. Therefore, one has the following 
perturbation equations from the first order terms 0{Sbi, Sp), 

3 3 

YSbi + 2YHB^6k = -^^kjSp-^^ktpBSp, (38) 



3 



Sk + ^Sk + HB^J2^^, = ^l-lklSp-^^k^^pB6p, z = l,2,3, (39) 

3 

Sp + 3-fHB6p + jY^^'PB = 0- (40) 



In order to solve the above system of differential equations, we need an inspiration from the process of constructing 
the exact solutions. First of all, from the general results (15, 1^) in both the GET and the brane theory, their dynamics 
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on the perturbation variables follow the constraints 

3 

5p = ^-ipB^Sh,, (41) 
1=1 

3 3 

5k = \ Y.^b,-^iZ^r (42) 



3^ . , 



By summing the Eqs. ( p39[ ) and then using the result (41), we end up with a second order differential equation for a 
variable of the combination 

3 3 3 

^ 5k + 6Hb 5k - 7^1 PB [3(7 - 2)fc2 + (7 - l)^^^^] ^ - 0. (43) 



The task is to solve ^ 5bi from the above equation and then to construct 5bi and 5p from the Eqs. ( [42| ) and (^ij). For 
the purpose of the stability analysis in the final stage, it is sufficient to consider the large time limit behaviors of the 
perturbation variables 5bi and 5p. From the discussion on the asymptotic behavior of the exact solutions, qualitative 
outcome in both the GET and the brane theory, we should divide our analysis into two different cases in the presence 
of absence of the cosmological constant. 

For the case with a positive cosmological constant, one can extract the asymptotic forms of the background variables 
from the exact solutions which give 

Hb^\/A/S, Vb (x e^piVsAt), ps cx exp(-\/3A7i). (44) 

The asymptotic expression of indicates that the third term in Eq. (|43| ) can be neglected in the large time limit. 
As a result, we 1 
the brane world 



As a result, we have the following equation for the asymptotic Xli=i which remains valid in both the GET and 



J2 + %/T2A Y5k = 0. (45) 



This in term leads to the final result 



Y5bi cx exp(-%/T2Ai). (46) 
i=i 

Therefore, from Eqs. (^ ^2|) we can obtain the asymptotic expressions of the following perturbation variables 

(56, cx exp(-\/T2At), (5p cx exp[-V3A (7 + 2)i]. (47) 

This indicates that that the background solutions are stable. 
Similarly, for the case with A = 0, we have 

2/7 



Hb ^ — Vb 



( ySpokj 7 \ ^ 4 2 
o M ' TuTZ^ * ' (48) 



37 " \^ 2 J . 3^2^ 

where the preceding coefficients become important in this case. Therefore, the Eq. (^3|), when the brane correction 
due to the quadratic energy density term signified by the parameter can be neglected in the large time limit, 
reduces to 



33 3 
Y^5k + -t-^Y.^k-'^^^^^t-^Y.^b^ = ^- (49) 

The solution of the above equation has the form 

cxt", (50) 



i=l ^ z=l ' i=l 



3 



i=l 
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where the exponent parameter is determined by 



4 K^/ — 2) 

a{a-l) + -a-^ ^ = 0. (51) 

7 7 



The explicit expression of a is 



^^t^3±Vi^EW±S^tEl^,^ (52) 

27 

which is always negative for all possible values of 7. The only exception is the case when 7 = 2 such that a = is a 
possible solution. In summary, the asymptotic behavior of perturbation variables Sbi and Sp is 

5h(xe, Sp(xt"~'^, (53) 

which shows that the background solutions are always stable even when the cosmological constant is vanishing. 



IV. ANISOTROPIC UNIVERSE WITH SCALAR FIELD 



In the previous section we have considered the evolution of the universe with a perfect fluid, obeying a barotropic 
equation of state, in the anisotropic background. We will study the stability problem of the anisotropic universe with 
a scalar field in this section. It is known that the scalar field with a potential can be thought of as a perfect fluid 
with the energy density and pressure given by = 0^/2 + U{cj)) and = 0^/2 — U{(p) respectively. Here C/(0) is the 
scalar field potential. The scalar field is expected to play a fundamental role in the evolution of the early Universe. 

In this section, we will only consider the case where the scalar field potential is a constant, U{(j)) = A = constant > 
acting as a cosmological constant. As discussed earlier, the related anisotropic cosmology in the GET is just a particular 
limit of the brane universe. Therefore, we will consider the brane universe in details for the moment. The GET can 
be recovered by taking a suitable limit. 

When we discuss the stability problem of the system, we will focuss on the system with a scalar potential that 
admits at least a local minimum 4>o such that U{4> = 4>o) = A. We will assume that 0o is the asymptotic solution to 
the field equation in the large time limit. We will able to show that the system remains stable in the large time limit 
even the evolutionary solution is only known when the scalar potential is a constant. 

A. Brane Cosmology 

For a Bianchi type I brane Universe with a scalar field the gravitational field equations take the following form : 



-kf 



^J^i^H,) = klU{cb)-^kU^-U^{cb)y z-1,2,3. 



(54) 
(55) 



Note that the scalar field also obeys the following evolution equation: 

+ 31/0+^=0. 



(56) 



By imposing the constraint U{(f>) = A, the Eq. ( |5^ ) gives = 2(I)q/V with 0o > a constant of integration. Summing 
the Eq. (|5^) over the index i and comparing with the Eq. (^4|) , one can extract the following equation: 

V ^ K^V - kIV-^. (57) 

Here we have denoted the positive constants k and kq as = 3fc|A + fc|A^/4 and Kq — fcgc^g. 
The general solution of the Eq. (^) can be shown to be [Q 



V{t) = i-e-«(*-*«) sj [e2«(t-*o) -Cf - AkIk^ 



(58) 
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Here C is an integration constant. 

The time evolution of the expansion, scale factors, mean anisotropy, shear and deceleration parameter are given by 



H{t) 
ai{t) 
A{t) 



3 



(e2«(t-to) _ cy 

a^^e-«(*-to)/3 

12i^2g2«:(t-to) 



\ 2 


1/6 






exp 


2K^F 



1,2,3, 



„2K(t-to) 



cy 



12e 



re4«(t-to) _ C2 + 4k2k21^ 



(59) 
(60) 
(61) 
(62) 



where aoi, i = 1, 2, 3 are arbitrary constants of integration and F{x) = /[(a^^ — C)^ — 4k§k2]-i/2(^2;. We do expect 
that the Universe starts its evolution from a singular initial condition. Therefore, according to the value 1^(0), the 
parameter tg should be chosen as exp(— 2k<o) = C + 2kok. 

Note that we can analyze the evolution of the Universe in two different regions. First of all, at the very early time, 
the parameters can be shown to be 



V ~ V2^ - 0, QKH/ko - 0, g ~ 5 > 0. 

Secondly, the asymptotic behavior of the parameters can be shown to be 

Vo^e^K Acxe-^f^* o ^ -1< 0. 



(63) 
(64) 



This result is similar to the model with a perfect fluid such that the brane Universe both evolves from an isotropic 
singular state initially. Later on, the mean anisotropy parameter A increases dynamically and decays to zero after 
the mean anisotropy parameter reaches its maximal value. 



B. Conventional Einstein's Theory 



The conventional Einstein's theory (GET) can be easily deduced from the above results by taking the limit = 0, 
or, in this case, kq — 0. The result is 



V{t) 
Hit) 
ai{t) 
Ait) 



c 

~2k 



K e 



2Kt 



1 



3 e2«* - 1 ' 



Kni/3. 



exp 

-2 



i = l,2,3. 



(65) 
(66) 
(67) 

(68) 
(69) 



Note that exp(— 2Kto) = C. The asymptotic behavior is the same as the case in the brane world. The initial state is, 
however, very different such that 



V . 



Ct - 0, 



A - iK'^/C constant, g 2 > 0. 



(70) 



Therefore, for the GET, the Universe has to start out from an anisotropic initial expansion and then decays all the 
way to the phase of isotropic expansion in the large time limit. 

The evolution of the mean anisotropy parameter A{t) for the GET and the brane world is presented in the Fig.|^. 
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FIG. 5: Mean anisotropy parameter for the Bianchi type I universe with a scalar field: for the GET (solid curve) and for the 
brane world (dotted curve). The normalization of the parameters is set as k = 1, 2koh = C. 
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FIG. 6: Deceleration parameter for the Bianchi type I universe with a scalar field: for the GET (solid curve) and for the brane 
world (dotted curve). The normalization of the parameters is set as k = 1, 2kok — C. 

C. Stability Analysis 

Since the the quadratic brane correction plays a significant role only at the early stage of the evolution of the 
Universe. For the purpose of our stability analysis in the large time limit, we will focuss on the asymptotic behavior 
of the fields. Therefore, it is sufficient to consider the background solutions in the conventional Einstein's theory. 

Perturbations of the fields of a gravitational system against the background evolutionary solution should be checked 
to ensure the stability of the exact or approximated background solution. In principle, the stability analysis should 
be performed against the perturbations of all possible fields in all possible manners subject to the field equations and 
boundary conditions of the system. In the following section, we will divide the perturbations into two disjoint classes: 

(a) the perturbations of the scale factors, or equivalently the metric field; and (b) the perturbations of the dilaton 
field. 

We will argue that the most complete stability conditions we are looking for can be obtained from class (a) and class 

(b) perturbations; even the backreaction of the scalar field perturbation on the metric field perturbations is known 
to be important We will show that this backreaction does not bring in any further restriction on the stability 
conditions. 

The reason is rather straightforward. One can write the linearized perturbation equation as 

i^l/a, + Z?;(50 = 0, (71) 

for the system we are interested. Moreover, perturbations are defined as ai — + dui and (/) = (po + S'P with the index 
denoting the background field solution. Note also that the operators Dl^. and denote the differential operator 
one obtained from the linearized perturbation equation with all fields evaluated at the background solutions. The 
exact form of these differential operators will be shown later in the following arguments. 

One is looking for stability conditions that the field parameters must obey in order to keep the evolutionary 
solution stable. One can show that class (a) and class (b) solutions are good enough to cover all domain of stability 
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conditions. Let us denote the domain of solutions to class (a), (b), and (a+b) stability conditions as S{a),S{b), 
and S{a + b) respectively. Specifically, the definition of these domains are defined by S{a) = {6ai\Dl^,Saj = 0}, 

S{b) = {d(l)\D'^S(j) = 0}, and S{a + b) = {iSai,d(j))\Di^6aj + DiS(j) = 0}. 

Therefore, one only needs to show that S{a) n S{b) C S{a + b). This is because that ^'D]^.6aj = and D^^Scf) = 0" 
imply that "Dl^Saj + D'^Scj) = 0." On the other hand, ''Da^Saj + D'^Scj) = 0" does not imply that "DI^Sqj = or 
D]jj5(t) = 0." Hence class (a) and class (b) solutions cover all the required stability conditions we are looking for. 
Hence we only need to consider these two separate cases for simplicity. 

In addition, one knows that any small time-dependent perturbation against the metric field is known to be equivalent 
to a gauge choice j2^. This can be clarified as follows. Indeed, one can show that any small coordinate change of the 
form x'^ — — will induce a gauge transformation on the metric field according to g'^^ = g^^ + D^e^ + D^e^. 
Therefore, a small metric perturbation against a background metric is amount to a gauge transformation of the form 
a[ = ai -\- e^'hi for the Bianchi type-I metric with = (constant, ei(i)). This is then equivalent to small metric 
perturbations. If a background solution is stable against small perturbation with respect to small field perturbations, 
one in fact did nothing but a field redefinition. 

If the background solution is, however, unstable against small perturbations, e.g., the small perturbation will 
grow exponentially as we will show momentarily, the resulting large perturbations cannot be classified as small gauge 
transformation any more. Therefore, the stability analysis performed in the literature ||l^, |l^, |2^, ^ ^ 
for various models against the unstable background solution served as a very simple method to check if the system 
supports a stable metric field background. This is the reason why we still perform a perturbation on the metric field 
for stability analysis; even a small perturbation is equivalent to a gauge redefinition. 

Note that one should also consider a more general perturbation with space perturbation included. The formulation 
is, however, much more complicated than the one we will show in this paper. We will focus on the time-dependent 
case for simplicity in this paper. The space-dependent perturbation analysis is still under investigation. The time- 
dependent analysis alone will, however, bring us much useful information for the stability conditions about the model 
we are interested. For example, we will show in the following subsection that the solution found in Ref. remains 
stable as long as the scalar field falls close to any local minimum of the potential U{(j>). Note again that the solution 
found in Ref. iQ is an exact solution only when U ~ constant. 

Dilaton Perturbation 

Let us consider the perturbation of the dilaton field of the following form 

(/) ^ -f (50, C/(0) ^ A + d^8<l>. (72) 



By setting fcs = in the Eqs. (M, 55, M), the field equations for become 



- 20B<50 + a^C/(5(/) = 0, (73) 
d4,UH = 0, (74) 
■iHB5^ + dlU5(l) = 0. (75) 



Due to the fact that Hbi -5 , the asymptotic solution of S(t) can be shown to be 



(5(/)cxexp {-K±y/K^ -4/3)i/2 . (76) 

Here /? = d'^U\t^oo- Therefore, the perturbation of dilaton field decays to zero if > 0. Hence we show that the 
system is stable with respect to the scalar perturbation if /? > such that the scalar potential is capable of confining 
the scalar field to its local minimum. 

Note that there are two additional equation [/' = and (j) — which is required for consistency of the stability of 
the system. These constraints simply imply that the scalar field must be at rest at the local minimum of the scalar 
field potential. As we have pointed out earlier that the background solution we have at hand is an exact solution for 
the model with a constant cosmological constant. Therefore, the background solution we used for stability analysis is 
only an approximated solution which remains valid only when the scalar field is close to the local minimum of the the 
scalar field potential. As a result, one should not be too serious about these two further constraints. In fact, these 
constraints are both negligible in the large time limit when the scalar field falls close to the local minimum of the 
scalar field potential, namely, ?7' — > and ^ ^ 0. Hence one is able to show that these additional constraints can be 
made satisfied approximately in the large time limit. 



Metric Perturbation 
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Using the metric perturbation (p7|) , the perturbation equations for the metric perturbation 5bi can be obtained 
from perturbing the Eqs. (jsj, 55, 5^. The result is 

J2^b^ + '^J2^B^Sk = 0, (77) 

i i 

Sb, + -^Sb, + HB^y 6bj = 0, (78) 
Vb 

J 

^sJ^Sk = 0. (79) 

i 

Here we also choose the Hmit — 0. Note that the background variables Vb, clbz cind Hsi approach 

Therefore, the asymptotic behavior of the metric perturbation Sbi can be found from the second equation. In addition, 
the other two field equations will provide a constraint for the system. The results are 

-> c,e-«*, ^c, = 0. (81) 

i 

Hence, the metric perturbation 

(5ai = aB»(56, 6^2"*/^ 0, (82) 

in the large time limit. This indicates that the background solution of the system is stable against the metric 
perturbation as shown above. 



V. CONCLUSION 



We have discussed the anisotropic property of cosmological models in the GET and brane theory. A realistic model, 
being consistent with the current observations, should produce a small value of anisotropic parameter at the later 
stage of the evolution of our Universe near the last scattering surface. By assuming the Bianchi type I space-time for 
the evolution of our Universe, we found that the final state of the evolving Universe always approaches the phase of 
isotropic expansion in both theories. 

These two different theories give completely different initial anisotropy at the very early stage of evolution. Indeed, 
for the GET, the anisotropy tends to be large in the very early stage. In another words, the universe tends to begin 
from a highly anisotropic initial state. The mean anisotropy parameter A will then decay to zero as the time increases. 
On the other hand, the early time behavior of the Universe in the brane world scenario changes significantly due to 
the quadratic correction on the brane. As a result, any non- vanishing mean anisotropy parameter, A(t), tends to 
vanish in the very early period. There is a characteristic time, tc, that divides the evolution of A(t) into two different 
stages. The mean anisotropy parameter is increasing when t < tc and reaches its maximal value at i = tc- After that, 
A{t) starts to decay. This kind of behavior is clearly shown in the Fig|| and Fig.^. This result remains true for both 
the model with a perfect fluid and the model with a scalar field. And this appears to be a general feature independent 
of the types of matter considered. 

It is worth noting that the only exception is the model with p = (i.e. 7 = 1) of the perfect fluid model. The mean 
anisotropy of this model behaves similar to the models in the GET where mean anisotropy parameter is large in the 
very early time. Moreover, we also analyzed the stability problem for those exactly solved anisotropic models shown 
in this paper. The result indicates that all of the solution known to us are stable in the large time limit. Therefore, 
the evolution of the Universe in the GET starts with highly initial anisotropic expansion. The dynamic of the system 
will take the Universe to the phase of isotropic expansion in the large time limit. We also show that the final isotropic 
expansion will remain stable in the large time limit. In addition, the mean anisotropy parameter will keep decreasing 
as time increases. The model provided here is a useful and explicit model that is capable of providing us with a 
Universe that has a tiny anisotropy left over near the last scattering surface. 
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